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Abstract 

Galileon models, which were developed in the context of modified gravity, give a class of La- 
grangians containing derivative interactions without introducing higher order derivatives in the 
equations of motion. Here we extend the analysis to an arbitrary number of scalars, and examine 
the restrictions imposed by an internal symmetry, focussing in particular on SU(N) and SO(N). 
This therefore extends the possible gradient terms that may be included when constructing topo- 
logical objects such as sigma model lumps. 



* antonio . padilla@nottingham. ac . uk 
^ paul. saffm@nottingham. ac . uk 
' ppxsyz@nottingham.ac.uk 



1 



I. INTRODUCTION 



In an attempt to model the generic features of a particular class of modified gravity 
scenarios , Nicolis et al [EJ recently developed the notion of the galileon. The galileon 
is a scalar field, 7T, whose dynamics is described by a Lagrangian that is invariant under 
Galilean transformations of the form 7r — > n + h^ 1 + c, where and c are constant. A 
celebrated example of a galileon theory is derived from the DGP model To see this 
one can compute the boundary effective theory and take a limit in which an additional 
scalar degree of freedom decouples from the graviton fluctuations Q. This scalar can be 
identified with the position of the DGP brane in the Minkowski bulk, and its dynamics is 
described by a Galilean invariant theory in the decoupling limit. The Galilean symmetry is a 
descendent of Poincare transformations in the bulk. Most of the interesting dynamics occurs 
in this scalar sector, and one can study a number of phenomenological issues including, self 
acceleration jgj], ghosts [§], and the Vainshtein mechanism (Io| . 

Taking the Galilean symmetry as their starting point Nicolis et al [E| showed that if we 
also require the field equations to be at most second order in derivatives, then a generic term 
in the single galileon Lagrangian is of the form 

rf\^--V^ Vm] *d^---djr m ir, (1) 

where m < d, the number of spacetime dimensions. Many aspects of the single galileon 
theory have been studied, ranging from phenomenology [llj], covariant completions [Hj], an d 
their interpretation in terms of the position of probe DBI branes jl3j. See [l5[ for early 
related work. 

The extension of this symmetry to two galileon fields has been recently developed [14| 
(see also [3]). As with the single galileon this was motivated by modified gravity scenarios, 
and in particular, co-dimension two braneworld mo dels H3-E3. The two galileon fields can 



be identified with the position of the 3-brane in a six dimensional bulk. Generically there is 
no reason to expect there to be an internal symmetry linking the two brane bending modes 
since the brane can fluctuate independently along the two orthogonal directions normal to the 
brane. However, one could speculate that the decoupling limit of a rotationally symmetric, 
co-dimension two brane would correspond to a bi-galileon theory with an internal SO (2) 
symmetry - the Galilean symmetry coming from bulk Poincare invariance, and the SO (2) 
from cylindrical symmetry that remains unbroken by the brane. 

Here we take a slightly different view from previous work on galileons in that we do 
not concern ourselves with gravitational physics. Instead we will consider galileons purely 
from a field theorist's perspective. The Lagrangians describing galileon theories typically 
include higher order gradient interactions. Higher order gradient terms play a central role in 
constructing sigma models with stable lumps, such as the Skyrme [22|, 23[ and baby Skyrme 



model [24(, as they behave differently to the canonical kinetic term under a rescaling of co- 
ordinates. This allows one to evade Derrick's theorem 25|, 26], and so construct finite energy 
topological defects. These solitons have been used to model baryons, so naturally it would be 
interesting to extend these models, perhaps by introducing new types of gradient interactions 
to stabilise new configurations, or at least render them static. However, generically, such 
interaction terms will introduce the ghost-like instabilities associated with higher derivative 
field theories [23]. The nice thing about galileon theories is that even though higher order 
gradient terms are present in the Lagrangian, the equations of motion are no more than 
second order in derivatives. Thus the ghost problem is avoided, and one is encouraged to 
extend the Skyrme model in new ways. 
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In this paper we keep the general structure of the galileon term, but we allow for the 
scalar to be a multi-component field. If we were not to impose any symmetry between 
the components this would lead to a huge proliferation of terms rather quickly, but by 
imposing an internal symmetry we are able to find a manageable set of terms. We start in 
sections [II Al and lH Bl by constructing terms that utilize real scalars with an SO(N) symmetry, 
before moving on to a complex scalar realizing SU(N) symmetry in sections III CI and III Dl 
In section IHH we discuss Derrick's theorem in the context of galileons. By means of an 
example, we show explicitly how Derrick's theorem can be evaded, and one can indeed use 
the higher order gradient interactions to construct static soliton solutions. This result is in 
contrast to 28|, where it was argued that galileon interactions could not be used in this way. 



The difference here is that we can consider multiple fields constrained to lie on a non-linear 
manifold, and so one cannot consider field deformations orthogonal to the pre-defined target 
space We will summarize our results in section IIVI 



II. MULTI-COMPONENT GALILEON 

A general multi-galileon theory, in four dimensions, with A" real degrees of freedom is 
given by the Lagrangian (TiJ 

5 

4r = E « H - im ^d^TT^.d^n^, (2) 

m=l 

where {a n '"* m } are free parameters of the theory and ^r^ '^j = m ^ v t x \ Vl --- v t m v y Note 
that throughout this paper summation over repeated Lorentz (Greek) and galileon indices 
(Latin) are understood. Note further that we define the m = 1 term of expression ($Z§ to 
be a^n^. For each galileon term (ie, terms grouped by a a ll " Am ), since {/i/c} and {uk} 
are anti-symmetrized respectively, d flk d Uk can be "moved" by integration by parts from 7ii k 
to 7Ti 15 thus the galileon indices of ^^l '^ m y^i 1 d^ 2 d V2 iri 2 ...d^ m d Um Tti m are symmetric, which in 

turn means that we can choose a !l,,,tm as symmetric. The total number free parameters in 
a general A^-galileon theory is therefore given by 

J^(N + m-l\ _* (N + m-l)\ 

^ V m J 2^ m HN- 1)! ' 1 j 

m=l x ' m=l v ' 

For example, when A^ = 3 there are already 55 free parameters. 

Now, if there is an internal rotational symmetry within the A^-galileon, we will see that 
the number of free parameters is greatly reduced. 

A. SO(N) fundamental representation 

Suppose that the A^-galileon transforms according to the fundamental representation of 
SO(N), that is, 

- >()/::,. 0+0^ = 5^ k,l — 1, N. (4) 
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In order for the terms in ([2]) to be invariant under such a transformation we require the 
parameters {a n ""* m } to be invariant tensors of the fundamental representation, 

n ,\ : ■'■ (), h .. A) ; im . (5) 

Invariant tensors can be constructed from the primitive invariants of the group in the 
corresponding representation; the primitive invariants of the fundamental representation 
of SO(N) are the Kronecker delta and the Levi-Civita symbol [29| 

5 lj , e h - iN . (6) 

Since a M """ lm is symmetric, we can ignore the anti-symmetric Levi-Civita symbol and con- 
struct a n - lm solely from 8 l K Then it is clear that, up to 5th order, the non-zero free 
parameters are and a^ kl , which can be parametrised as 

a ij = a {2) 5 ij , a ijkl = a {i) 5 {ij 5 kl) . (7) 

So the most general SO(N) fundamental galileon Lagrangian is given by 

= a m vr J nvr, + a (4) tf£f i d p & r >K i d ll &'>K i d x d r >K j (8) 
= a m t^Ott + a (4) ^ {fd^Tx) {d,d^d x d T v) . (9) 



B. SO(N) adjoint representation 

We have seen that constructing a general SO(N) galileon Lagrangian for any N ^ 2 in the 
fundamental representation is straightforward, but it is not necessarily a simple procedure 
for scalars that transform in representations other than the fundamental. Another possibility 
we are particularly interested in is the adjoint representation, which may be conveniently 
described by iV x N anti-symmetric matrices, LT^ = — ILjj, transforming as 

n^ono T , or Uij^o^o/, (10) 

where O = (O i J )NxN again is a SO(N) defining matrix. Therefore a general SO(N) La- 
grangian with a galileon transforming in the adjoint representation may be written as 

5 

Cti = J2 P hjl '- ,imjm V [ : ] ^nnddU i2h ...ddU imjm . (11) 

m=l 

Note that for simplicity we suppress the Lorentz indices in this section. Treating the pair 
{ikjk) as "one index", ^ l «i>-,w'm can De considered as adjoint invariant tensors. Since 
ri^"^ddlli 2 j 2 ...ddlli m j m is symmetric in the galileon pairs, p n ^^--^ m can then be chosen as 
symmetric in the galileon index pairs (within each pair, exchanging i^ and jj~ is of course anti- 
symmetric). To determine all possible / 9 l i-7'i>--->Wm ) we would like to find all the symmetric 
primitive symmetric invariants up to 5th order (in terms of pairs) in this representation, and 
then construct ^u'li-iWm us j n g these primitive invariants. As an aside, since the indices (not 
pairs) still transform under fundamental SO(N) matrices, we expect, as we shall see shortly, 
that all the primitive invariants in the adjoint representation can actually be constructed 
from and e ll ' ijv . 
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Before the construction, however, we would like to clear up some special cases of the low 
rank groups. 5*0(2) is Abelian, and not surprisingly the adjoint representation of S0(2) 
has only one independent component, thus the adjoint SO (2) galileon reduces to the single 
galileon case. The adjoint of SO (3) is equivalent to the fundamental, so the adjoint of 
SO (3) offers no new model. 50(4) is not a simple Lie group and is locally isomorphic to 
SU{2) x SU(2), so the adjoint of 50(4) is reducible. Nevertheless, in terms of constructing 
Lagrangian singlets, we shall see that it goes along similar lines as some other simple 50(2/). 
We shall ignore 50(2) and 50(3) and consider SO(N) with N ^ 4 in the following. 

For a simple Lie group G, it is well known that there are /(= rank(G)) symmetric primitive 
invariants in the adjoint [30] . For simple 50(2/ + 1) and 50(2/), / and / — 1 of them are given 
by the first / and / — 1 even order (in terms of pairs) of the /c-family symmetric invariant 
tensors respectively, defined by the symmetrized trace 

fJljl,---imjm — \ STrfT 1 ^ 1 ^ 1 jri-mjm} (12) 

ml 

where T %k ^ h are the generators of the fundamental representation of SO(N), conventionally 
chosen as 

{T kl ) i 3 = -i (5 k 5 l1 - 5\5 kj ) . (13) 

By noting that the IT terms in the Lagrangian (ITT1) are already symmetric under interchang- 
ing the pairs {ikjk), the symmetrizing of the trace is done for us, leading to the most general 
Lagrangian for these cases to be 

£ n = /3 {2) Tr(nnn) + /3 ( 4)77 [ ; ] Tr(naan)Tr(aanaan) + ^rjQTriUddUddUddU). (14) 

For 50(2/) there is one more primitive invariant, associated to the Pfaffian, which for 
/ = 2,3,4,5 allows us to construct terms that may appear in the Lagrangian; for larger / 
they lead to terms with more than five lis, and so cannot contribute. The invariants are 

1 = 2: e ujii2i 2; ( 15 ) 

/ = 3 : e nji-w 3) ( 16 ) 

1 = 4: e hh-i4j^ ( 17 ) 

1 = 5: £ iiii-ius ) (18) 

which we note are symmetric in moving pairs of indices. The existence of these invariants 
is related to the fact that the determinant of an anti-symmetric matrix (Mij)2Nx2N can be 
expressed as the square of the Pfaffian, 

Pf(M) = y/det{M) 

= Jih-iNJN /if. . jjf. . (iq) 

2N]\f\ lvl i\n--- lvl iN3N- \ LJ ) 

These invariants give rise to additional terms for each of the four groups (50(2/), / = 
2, 3, 4, 5) respectively: 

£4 = m e(nnn) + e(n<9<9n)Tr(<9dn<9<9n), (20) 

£6 = ^ -» e (UddUddU) + ifffl^ e(UddUddU)Tr(ddUddU) : (21) 
£8 = mv - e (UddUddUddU) } (22) 
£io = riior]^ eiUddllddllddUddU), (23) 
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with the definition 



e(UddU...ddU) = e hj ^ jm U iljl ddU i2h ...ddU imjm . (24) 
Note that for 5*0(4) one might have expected an additional singlet term proportional to 

i]Qe(IlddIl) e(ddUddU). (25) 

However, by noting the relation 

we see this term degenerates to a linear combination of f7["']Tr (UddU) Tr (ddUddU) and 
^■^(UddUddUddU). 



C. SU(N) fundamental representation 



In this section, we promote the N galileons to be N complex fields fa and impose SU(N) 
symmetry within the complex iV-galileon = (fa, ...,<f>jf). The tensor product of two irre- 
ducible representations can produce a singlet only if they are conjugate to each other. For 
SO(N), the complex conjugate of the fundamental is equivalent to the fundamental itself, 
thus we could build singlets merely from the fundamental representation. For SU(N), this 
is generally not the case, so we include fields that transform with the complex conjugate 
representation 0* in building the SU(N) fundamental galileon. Since the conjugate repre- 
sentation of SU(N) is generally inequivalent to the fundamental representation, we would 
have to differentiate the indices of the two representations carefully. A simple way to do so 
is representing the fundamental field with a lower index fa and representing the conjugate 
field with a upper index fa = (fa)*, such that they transform as 

fa -> U % j fa, fa -> U'jfa, r, j U k j = 5t (27) 

A general SU(N) fundamental Lagrangian would be given by 

5 

= X) ah ''' lm ^ 2 2 :Z] faA^ 2 fa 2 -d„ m d Um fa m + terms with conjugate reps, (28) 

m=l 

where "terms with conjugate reps" represents terms with all possible combinations of the 
fundamental and conjugate representations, ie, all possible terms with upper in (corre- 
spondingly lower Zfc in a'" ). The primitive invariants of SU(N) fundamental and conjugate 
representations are given by the Kronecker delta and the Levi-Civita symbols j29| 

<5}> e 1 N 7 e ii...ijv! (29) 

and, as with the SO(N) case, we can construct our invariant tensors or" out of them. Since 
a " are totally symmetric, i.e., 

U JU2- " 3132- 31 32- \ OU ) 
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we can ignore the anti-symmetric Levi-Civita symbols and construct or" solely from Sj. 
Therefore, up to 5th order, the only non-zero free parameters are 

a y = a (1A) S], a ik jl = a {2i3 /$\ (31) 
and then the SU(N) fundamental galileon Lagrangian is given by 

= «(!,!) pOfri + a m v^d^&d.d^d^^ (32) 
= a (1)1) 0t D + a{2 2) ^ (^d p d°<P) (dpVtfdxdrtf) . (33) 

D. SU(N) Adjoint Representation 

We now turn to construct the most general Lagrangian for a galileon transforming under 
the SU(N) adjoint representation. Similar to that of the SO(N) adjoint, it is convenient 
to use A x A traceless hermitian matrices $ = ($/)atxat to represent the adjoint SU(N) 
galileon, which transforms under a SU(N) adjoint action as 

$ U$rf (34) 

where U = (Z7/)jv x jv again is a SU(N) defining matrix. Therefore a general SU(N) adjoint 
Lagrangian may be written as 

£* = E ^ V^rfi K hdd ^ h - dd ^J m - ( 35 ) 

m=l 

Again we suppress the Lorentz indices in this section. Treating the pair (ikjk) as "one index" , 
(3 n • "' lm ■ can be considered as symmetric invariant tensors (within the pairs, i& and are 



hermitian). So the task again is to find all the symmetric primitive invariants up to 5th 
order (in terms of pairs) in this representation and then construct /3 H j 1 " lm j m from these 

primitive invariants. Since fi 11 ^ " lm ^ always have upper and lower indices and products of 
e n "- tN and £j 1 ...j N can be reduced to product of Sj, we would expect that all these primitive 
invariants can be solely constructed from 5j. 

SU(N) is a rank N — 1 simple group, so there are N — 1 symmetric primitive invariants, 
which may also take to be given by the symmetric traces of the adjoint generators 



k\z! m jm m = 2,...,N. (36) 



The adjoint generators of SU(N), however, are more complicated than those of SO(N), 
including N(N — l)/2 the anti-symmetric generators, essentially the same as that of SO(N), 
as well as (N + 2) (A — l)/2 symmetric ones. Nevertheless, since the only building block 
we can use is Sj, it is not difficult to construct a new family of invariant tensors suitable for 
constructing the B n - "' %m ■ : 

V\r lm , =5^...5Y ) . (37) 

The first A — 1 (m = 2, A) of P %1 j u '".* m j m are primitive. Indeed, P %1 j u '".* m j m gives rise to 
the m-th order trace r]< "\Tr($>...dd$>). Therefore, for SU(N) with A ^ 5, the most general 
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Lagrangian with galileon fields transforming under the adjoint of SU(N) is 



£ $ = 72 Tr($n$) + y 3 riQTr($dd$dd&) 

+ y 4 riQ r Ti($dd$dd$dd$) + y A r)£ ] Tr($dd$)Tr(dd$dd$) 
+ ~f B 7iQTr($dd$dd$dd$dd$) + ^■■ ] T±($dd$)Tr(dd$dd$dd$). (38) 
For lower rank groups we have degeneracies, as a result of the following identity for an 

■ ■ ■ 

h— 3m 

5] 1 ...5! m 1 A il , - im , =0 for m > N, (39) 



arbitrary tensor A 11 ^'" '"' 



where the group indices run from 1 to N. Using this result, we can derive the following 
identities: 



for N = 2,3, 4. 



Jd 4 x 7/ [ ;" ] T!r($(95$) 4 ) = Jd 4 x -7/ [ - ] r &($55$) r &((55$) 3 ) (40) 

y*d 4 x 77 [ - ] Tr($(aa$) 3 ) = y"rf 4 x ^ [ ;;; ] Tr($aa$)Tr((aa$) 2 ) (41) 
for iV = 2, 3 only. 

Jd 4 x rj^Tr^idd®) 2 ) = (42) 

for N = 2 only. 

Therefore, for S , f/(4), the 75 term degenerates with the 75 term, thus can be dropped. 
For SU(3), we can further drop the 74 term. For SU(2), we can further drop the remaining 
odd order terms, left with only the 72 and 74 terms. This is what one may expect, since 
SU{2) is locally isomorphic to 50(3). 



III. DERRICK'S THEOREM AND EXAMPLES 

In carrying out this work we have largely been motivated by the possibility of using 
these terms to find static top olog ical defects, and one of the key ideas to understanding 



this is Derrick's theorem [25|, |26(. Derrick's theorem uses the fact that static solutions 
extremize the energy functional, and so by considering rather simple deformations of a 
putative static solution, one is able to make statements about their (non) existence. An 



extension of Derrick's theorem to the single galileon case was considered in [28J, where they 
concluded that such terms cannot stabilize solitons, here we see how one can get around 
this conclusion. The key to understanding the statement in [28] is that they were thinking 



of the galileon as a linear sigma model and so considered the deformed field configuration 



</>\,u(x) = ^0o(Ax), 



(43) 



where A and u are constants, and 0o is the putative soliton solution. It is then a straight- 
forward calculation to show that the term with n scalars contributes the following to the 
energy functional, 



E n (\,tu) 



X d+2 



WE®. 



(44) 
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It is now clear that the energy functional is not stationary under variations of u. Indeed, 
if we reduce u all the terms in the energy functional will decrease. However, in non-linear 
sigma models, such as the Skyrme and baby-Skyrme models, we force the scalar field to 
live on some particular manifold, S 3 and S 2 in these respective cases. This means that we 
are simply not allowed to do the full rescaling put forward in (143 p . and we are restricted to 
u — 1, giving an energy for the rescaled solution (in four spacetime dimensions) of 

E(X) = X-'E^ + XE® + X 3 Ef + A 5 £f . (45) 

This has a minimum as we vary A, so solitons in non-linear sigma models are not ruled 
out. Note that these statements do not rely on any particular representation, or symmetry 
group. In what follows we give an explicit example of a static soliton using the fundamental 
representation of SO (4). 

The SO (4) non-linear sigma model Lagrangian is given by 

C = -i^vr^Tr - V[ Z X T] (vr r <9 p <9 CT vr) (^vr^vr) , (46) 

where we have made a choice of units such that the dimensionful parameter in front of the 
quartic piece is absorbed into a co-ordinate scaling. The corresponding Hamiltonian for this 
theory is given by 

H — j d 3 X ^± T ± - 25 { f £] [{^^{dadVd^TL) + 2^ Tl) {t T d b d e Tt 

+ \dal?d a lL + 8 { t C f\ (^dad'n) (d^d^n) , (47) 

where a, b, c, ... refer to spatial indices, and " dot" corresponds to differentiation with 
respect to time. We now define our static ansatz for the four-component field 

7r T = ( — sin fir), — sin f(r), - sin f(r), cos/(r) ) . (48) 

We will be interested in solutions with integer winding number n, satisfying f(r) — > 7m, as 
r — > oo. Substituting this into the Hamiltonian yields the following energy functional 

E[f] = 4vr fdr 2(f) 4 - 8sin/C ° S/ (/T + r4 + 24sin2/ 2 ~ 16sin4/ (/) 2 
J r 2r 

r 4 sin 2 f — 4 sin 4 f . . 

+ -> ( 49 ) 

where we have performed some integration by parts to eliminate f"{r). Varying the energy 
functional to get static solutions, SE/5f = , we obtain the smooth profiles for winding 
numbers n = 1 and n = 2 respectively, as shown in Fig. [TJ The energy of these two profiles 
is finite and, by explicit integration, take the values E^ ~ 101.8 and E^ ~ 462.3 in these 
units. Note that since E^ > 2E <yl \ the profile with winding number two is unstable and 
can decay to two solitons of winding number one. 

Finding static solutions with a topological winding number is one thing, but checking their 
stability is another. For example, although topological textures have a winding number they 
are certainly not stable, see for example [Hj]. So, while the above solutions are extrema of 
the Hamiltonian, they are not necessarily minima, and a full perturbative analysis of this 
system is left for future work. 
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FIG. 1: Static solitonic solutions for winding number n = 1 (solid-blue line) and n = 2 (dashed-red 
line). The spherically symmetic ansatz (|48j) is used. The energy of them are 

~ 101.8 and 

E(2) ~ 462.3 respectively. 



IV. SUMMARY 

In this paper we have expanded on the bi-galileon field theory studied in [13], which 
itself was a generalization of the galileon field theory examined in |5j . Motivated by certain 
braneworld setups and also to avoid rapid proliferation of possible multi-galileon terms, we 
have imposed internal symmetries for the multi-galileon field. There are clearly many differ- 
ent possibilities for the internal transformation properties. In this paper we have focused on 
four cases, namely the fundamental and adjoint representations of SO(N) and SU(N). We 
have written down the most general galileon theories with these symmetries, and shown that 
they contain a small number of terms. Incidentally, one can check the results using group 
theory arguments that we have found the correct number of terms 32| , the point being that 



these terms correspond to singlets in the symmetric products of the representations. For 
example, the adjoint representation of SU(5) is the 24; so, if we want to know how many 
($) 4 terms there should be, we need to consider 

24 x 24 x 24 x 24 = 1 © 1 © 24 © 24 © 24 © (50) 

ie, there are two singlets. This agrees with the number of quartic $ terms in (j38p . namely 
two. We have only explicitly shown results for symmetric multi-galileon in four spacetime 
dimensions, but it is clear how our techniques can be extended to any specific number of 
dimensions. 

In galileon models, as well as the canonical kinetic term, one finds terms with more than 
two derivatives. Such terms behave differently under a rescaling of coordinates, x^ — > Xx 11 , 
meaning that Derrick's theorem allows for the possibility of stationary finite energy solitons. 
Higher order terms in the action are not unusual in the theory of solitons, indeed, they are 
required in the Skyrme model. Here we are modifying the Skyrme model, constructing static 
solitons using these galileon terms. 

To have an explicit example, we have constructed hedgehog solutions for the 5*0(4) 
nonlinear sigma multi-galileon confined on S 3 . The static solutions we obtained possess 
topological charges, an encouraging sign for the stability of a soliton, but certainly not a 
proof, so we would like to conclude this paper by saying a word or two about stability. For 
a non-linear sigma model constructed out of a single galileon, it has been argued that zero 
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mode excitations inevitably suffer from ghost-like instabilities 28] . Those arguments do not 



necessarily apply here for the following reason. The act of constraining 7r to lie on an S 3 is 
equivalent to adding a constraint term A ( | tt | 2 — 1) to the Lagrangian. Perturbatively this 
gives rise to a "mass term" A(<57r) 2 , where the background Lagrange multiplier A depends 
on the 7r profile. Numerically, we can confirm that A is always negative; thus we come back 
to the same structure as kinks and skyrmions, where the zero-mode instability arguments 
of [28| do not apply. Of course, to confirm the stability, extensive numerical studies are 
required that are beyond the scope of this paper, and which constitutes our future work. 
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